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We present an analytical proof and a generalization of the Fayache-Sharon- 
Zamick relation between single particle energy splittings and the SU(3) limit 
q ' in Elliott's model. 
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I. INTRODUCTION 



In a recent publication p], Fayache, Sharon and Zamick compared the collective model 
result for a rotational band to the Od-ls splitting obtained with the OXBASH program ||. 
For the shell model calculations done with OXBASH they use a Hamiltonian consisting of 
a spherical harmonic oscillator and a residual interaction which is the isoscalar quadrupole- 
quadrupole (Q-Q) in coordinate space only, i.e., they omit the Q-Q term in momentum space 
and allow for interactions between shells with A A = 2. From these shell model calculations 
the authors of Ref. [jl| found that one third of the Od-ls single particle (s.p.) splitting (6%) 
comes from the interaction of the valence particle with the core while two thirds (12%) come 
from the diagonal Q ■ Q interaction. The Od-ls s.p. energy-splitting (18%) is the same as 
the energy splitting between L=2 and L=0 rotational bands obtained with Elliott's SU(3) 
model in the s-d shell @. A similar numerical relation was found in the f-p shell JE]. This 
is an interesting observation that deserves further study. 

In this paper we give an analytical proof of those numerical relations and show that 
they are particular cases of a general property of the Hamiltonian used. To this end we 
first derive complete analytical expressions for the s.p. energies (these will be defined more 
precisely in Sec. 2), and then derive the s.p. energy splittings in a major shell using the 
same Hamiltonian as in Ref. [1|]. We consider the problem of how to obtain the s.p. energy 
splittings in a major shell with such an interaction which will preserve the SU(3) results 

More specifically, we use the Hamiltonian 

ff= 2 ^£fc+*0-f£G(i)-(<w) (i) 

i ij 

where 

Q(i)-QU) = EQ^)QlU) (2) 

Q, (<) = b 2 x\ Kf fo) (3) 

with dimensionless coordinate and momenta Xj = r^/fr, 7Tj = Pib/h, and the harmonic oscil- 
lator (H.O.) length scale b = ^Jh/mu . 
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II. DEFINITION OF SINGLE PARTICLE ENERGIES AND ENERGY 

SPLITTINGS 



To parallel the shell model calculations of Ref. JT], we consider the effect of the residual 
Q-Q interaction when we have a closed core and put a valence nucleon in different orbital 
(N£m) states of a major shell (N) out of the closed core. In other words, we look for the 
s.p. energy splitting between Nl and N£' levels caused by the Q-Q interaction. To this end 
we write the s.p. energy of the Nl level as 

E m = hu ( N + I) + E%\ + E% (4) 
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where E®\ comes from the diagonal part of the Q-Q interaction (Q (i) ■ Q {%)) 

E% = -f < N£\Q ■ Q\N£ > (5) 



or 
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Ef e = -Ax<N£\x 4 \N£>, X = X^ (6) 

while E% comes from the interaction of the valence particle with the core. The direct term 
of the 2-body interaction (Q(i).Q(j),i ^ j) is proportional to the quadrupole moment of the 
core, hence it is zero for a closed core. Consequently E m is given by the exchange term 

Ef e = X E \<N£m \Q^\N c £ c m c >\ 2 (7) 

N c t c m c n 

Clearly the s.p. energy defined in Eq. (|4]) corresponds to the energy difference of the 
systems with A + 1 and A (core) nucleons, calculated as the expectation values of the Hamil- 
tonian ([!]) in the ground state of the unperturbed (x = 0) harmonic oscillator Hamiltonian, 
i.e., in the simplest, lower order harmonic oscillator shell model wave function. We adopt 
this definition to follow the numerical work described in Ref. although this is not the 
most general definition of s.p. energy that may be considered with the Hamiltonian ([!]). 

Using the explicit expression for H.O. wave functions in terms of Laguerre Polynomials 
we obtain the following result for E®] 



E%\ = -Ax {n{n - 1) + An (n + £ + - 
n + £ + £) (n + £+^ 
X 2£(£+l) -3(2N 2 + 6N + 5 



(8) 



with n the principal quantum number (n = (N — £)/2)(see for instance [[§]). 

To evaluate Eq. (0) we have to keep in mind that according to the selection rules for 
H.O. wave functions only matrix elements with N c = N — 2 will contribute to the sum 
over N c . This is clear since the Q M operator only connects the states with N to those with 
N' = N, N ± 2, and of these possible N' values only N' = N c = N — 2 belongs to the set of 
core levels. Carrying out the sum over m c and fi in Eq. (^]) we get 

2 / P 9 P \ 

Eft = 8xE (< N 1 \ x2 \ N - 2 4 >) (24 + 1) 







For H.O. wave functions the matrix element in Eq. (PI) is given by 



< N £ \x 2 \N -2£' >= 5 e/ j\jn(n + £+- ) +5 e4+2 \h> (n - 1! 



1 \ \n-£-^) (n + t-- 



(9) 



(10) 



Substitution of Eq. (|10D into Eq. (H), together with the following equalities that hold 
for 3 — j symbols, 



£(2f+l) 



( 1 2 e ^ 



(2£ + 3) 





^0 




^£ 2 £^ 


2 


f £ 2 £-2^ 




+ (2£ - 3) 




v y 




^00 



2£ + 1 



(11) 



(12) 



(2£- l)(2£-3) 



/ \ 2 
" 2 £-2 1 3£(£-r 
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2 (2£ + 1) 



(13) 



leads to the following result for the valence core interaction 

Ef e = 8x \n{n - 1) + n £ + %£{£ - 1) 



= x[^ + l) + 2iV(iV-2)] (14) 

It is now a simple matter to compute the energy splitting A N ^e,e') between iV-shell orbital 
partners N£ and Nt'\ 

A mn = {E m - Em') = SxW + 1) - ?W + 1)} (15) 

or for £' = £ - 2 

A w= <_ 2 ) = 6x(2£ - 1) (16) 

which agree with the SU(3) result. In particular in the N = 2 and N = 3 shells we recover 
the numerical results of Ref . (T| . 

We can also show easily that one third of the energy splitting A jv comes from the 

valence-core interaction and two thirds from the diagonal term of Q-Q. Writing 

Ajv(^') = A N \t t i>) + ^N 2 (i,P) ( 17 ) 

with 

AQl __ TTlQl 771 Ql / 1 Q\ 

^N(e,e>) = &m - &Ni' { 18 ) 

^N(e,e>) = ni - a Ni> K 1J ) 
we obtain using Eqs. (|S|) and (|I41) 

AS v) =2x [£(£ + l)-£'(£' + l)} (20) 

*%*)=TW + 1)-W+1)] (21) 

or 
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A%\ ijtf) = lA N(m I = 4(2£ - l)x for e = £ - 2 
Ag^ = lA mn | = 2(2£ - l)x for f = I - 2 

which again particularize to the numerical results found in Ref. for the N = 2 (ls-Od) 
and the iV = 3 (lp-Of) shells 

A% fi) = I2x A% A) = 20x 

Af 2fi) = 6x A? ( 2 3il) = lOx (23) 
A 2 ( 2 ,o) = 18x A 3(3,i) = 30x 
We note that although not explicitly mentioned, we have been considering the valence 
particle in a shell close to the core. In principle we may as well consider a valence particle 
in a higher iV-shell, > A^™ ax + 2. In that case the valence-core interaction is zero and 
only the Ql splitting remains, i.e., for these higher iV-shells one gets a smaller splitting. 

III. COMPARISON WITH ELLIOTT'S Q Q INTERACTION 

One may wonder how the one-third, two-thirds division of the s.p. energy splitting 
between core-particle and diagonal contributions might be related to the inclusion of the 
momentum-dependent parts of the Elliott quadrupole operator 



(22) 
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It is therefore instructive to study how the s.p. energy splitting is shared by the position, 
momentum, and crossed terms of the Elliott Q ■ Q interaction 

H QQ = ~^ EQ E d)-Q E (j) (25) 

/ a 

We stress that this interaction is formally identical to that in Eq. ([[]). The only difference 
comes from the replacement of the position quadrupole operator Q by the Q E operator. The 
latter is a sum of the dimensionless position and momentum quadrupole operators 

q^Mx^ix,) (26) 



# = yf ^im) (27) 

We recall that an important property of Elliott's quadrupole operator is that it has zero 
matrix elements between different N— shell states because the AiV = 2 matrix elements of 
q p exactly cancel those of q r . As can be seen from Table 1 the reduced matrix elements of 
are 



N't! 



Q E NZ) = 5 NN , x \^-b 2 \-8, e (N+l 



_ _ 1(1 + 1) (21 + 1) 

4tt" I V" ' 2)\ (2£-l)(2£ + 3) 



/3 (N - £ + 2) (N + £ + 1) £ (£ - 1) 
2 2£- 1 

^/f EEEgM±M±2 > (28) 

Therefore, with Elliott's interaction there is no particle-core interaction and the s.p. energy 
defined in the previous section [Eq. (Q)], only gets a contribution from the diagonal term 
defined in Eq. (||) with Q replaced by Q E . This diagonal contribution can be easily calculated 
from Eq. ( p8| ) and it is found to be 

Em = X [M {& + 1) - 4iV (N + 3)] (29) 

Hence, when Q is replaced by Q E in Eq. (U), the sum of the one-body (-Ejv<?) arLC ^ two-body 
fi^f J contributions to the s.p. energy [Eq. (§)] is replaced by the one-body contribution 
E m . The resulting value of the s.p. energy differs only in the N— dependence and therefore 
one is left with identical s.p. energy splittings: 

Kw) = E% E t - E%§ = 3 X {£ (£ + 1)- £' (£' + 1)] = A mn (30) 
We also note that this s.p. energy splitting can be decomposed in three contributions, 

one coming from the q r ■ q r interaction \A Er> j , one coming from the q p ■ q p interaction (^A Ep ^j , 



and one coming from the crossed q r ■ q p + q p ■ q r interaction yA Erp j. Using Table 1, it is 
straightforward to check that 



rpE r 
N£ 



-x(N£\q r -q r \m) = x 



)-£{£+!)-- (2N 2 + 6iV + 5) 



(32) 



EZ P e = -x(M\q p -q p \M) = E : 



Nl 



(33) 



E E N J = -X (N£ \q r ■q P + q P -q r \N£) = X 



2£ (£ + 1) - [N 2 + 3N - 



t) 



(34) 



Therefore the SU(3) s.p. energy splitting 3% [£ (£ + 1) — £' (£' + 1)] is shared as one-sixth, 
one-sixth, and two-thirds by the position, momentum, and crossed terms, respectively. 



A E r 



E * 



IV. A NOTE ON SPIN-ORBIT INTERACTION 



Since spin-orbit coupling is very important in nuclei it is also useful to examine how the 
above results are affected when the one body spin-orbit interaction (Vso = —^Xso £- s ), is 
added to the Hamiltonian. Different cases have to be distinguished. The simplest cases are 
the s-d shell and f-p shell nuclei where the closed cores consist of closed iV-shells (N = Z = 8 
or N = Z = 20). In these cases the results given in the previous section still hold for N£ and 
N£' orbitals. The only difference is that one may in addition consider the energy splitting 
between different £j and £'j' subshells. Since the energy splitting between (j — £ — |) and 
(j — £ + |) partners is not perturbed by the Q-Q interaction, the total N £ j- N £' f 
splittings are for £' = £ — 2: 



A N(£j,e>j>) = §x{2£ - 1) + Xso < 



=F2 for j — £ ± I and j' = j — 2 
-(2£-l) forj'=j-3 (j = £+k) \ ( 35 ) 
+(2£-l) for/=j-I (j = £-i) 
i.e., the spin-orbit in general destroys the SU(3) limit result, but for some of the £]-£']' 
splittings the proportionality to (2£ — 1) remains. 



On the other hand when the closed core corresponds to magic numbers with N or Z = 
28, 50, etc., where only one of the N£j subshells of the valence shell iV (the £j subshell 
with n = 0,£ = N, j = £ + |) is closed, one has to take into account in addition the Q-Q 
interaction of the valence particle with this subshell, i.e. when N or Z = 28, 50, 82, 126 the 
2-body interaction contains also the term 

SEf ej = X E l< N£jm\Q tl \N£ cJc m ]C >| 2 (36) 

rrij c fi 

with 4 = JV and j c = £ c + \ = N + \. 

This contribution gives an extra term to the energy of the valence subshells with quantum 
numbers 

(a) £ = 4 = JV, j = j c - 1 = N - \ (37) 



and 



(6) £ = £ c - 2 = N - 2; j = £ c - 2 + i = N - ~ (38) 



the value of this extra term is 



?Q 2 ^(2JV + 3)(JV + 1) 



(a) 8E%% = 6x ^ N - ^ m _ - 3x for N large (39) 



(b) 5E% = 24x iV 2 ( ^+ ] 1) - 12xiV for N large ( 10) 



Thus, this extra contribution may also spoil the SU(3) limit [Eq.(|T^)] for some of the 
s.p. energy splittings. 

V. CONCLUSIONS 

In summary, quite generally, the s.p. energy splitting in a major shell N preserving the 
SU(3) result 

&>N{t,e) = Eni — Ejsri' = Ajyj^ + A^ e ^ = Affy^ = 

= 3x \< N£\L 2 \N£ > - < JV£'|L 2 |JVf >1 (41) 



has been obtained with a Q-Q interaction in coordinate space that allows for interactions 
between AiV = 2 shells. As a general rule, the valence-core interaction — proportional to 
the matrix elements of the Q-operator between AiV = 2 shells — generates one third of the 
energy splitting in Eq.([y]), while the other two thirds are due to the diagonal (one-body) 
part of Q ■ Q. This result is surprising since traditionally the SU(3) limit of Elliott's model 
was derived restricting the action of the Q-Q interaction to a single major shell (eliminating 
the AiV = 2 matrix elements of the Q-operator). 

We have also explicitly shown that with Elliott quadrupole operator the same s.p. energy 
splitting is obtained by adding up the contributions from the position, momentum, and 
crossed position-momentum terms of the Elliott quadrupole-quadrupole interaction. This 
comparison is important because it illustrates how in this instance, the same effect can be 
obtained by taking into account (2-body) particle-hole interactions or (1-body) momentum 
dependent interactions. The fact that the s.p. energy splittings are equal while the s.p. 
energies are not, also serves to illustrate how the two types of interactions can be equivalent 
in some aspects while differing in others. 

Our conclusion here generalizes and reinforces the observation made in Ref. for the 
iV = 2, 3 shells and suggests new ways for further applications of the model Hamiltonian 
(H]). Since applications of this model Hamiltonian in the past have proven to be very fruitful 
it is worth to explore it in a new direction. A clear practical application involves deriving 
the Q ■ Q interaction from the realistic interaction. One now does not have to justify the 
momentum terms. 

We have also found that although the spin-orbit coupling tends to spoil the SU(3) limit 



in Eq. (f4TD, the proportionality to (2£ + 3) of the s.p. energy splitting between (£,j) and 
(£ + 2, j') subshells is still maintained in some cases. 



10 



ACKNOWLEDGMENTS 

This work was supported by DGICYT (Spain) under contract number PB95/0123 and 
by a Department of Energy grant DE-FG05-95ER-40940. 



11 



REFERENCES 



[1] M.S. Fayache, Y. Sharon and L. Zamick, Phys. Rev. C 55, 1575 (1997). 

[2] B.A. Brown, A. Etchegoyen and W. D. M. Rae, The Computer Code OXBASH, MSU- 
NSCL report number 524, 1992. 

[3] J. P. Elliott, Proc. Roy. Soc. London A 245, 128 (1958); 245, 562 (1958). L.S. Gradshteyn 
and I.M. Ryzhik, Table of integrals, series and products (Academic Press inc., London), 
1980; P. Morse and H. Feshbach, Methods of Theoretical Physics (McGraw Hill Book 
Company), 1953. 

[4] IS. Gradshteyn and I.M. Ryzhik, Table of integrals, series and products (Academic Press 
inc., London), 1980; P. Morse and H. Feshbach, Methods of Theoretical Physics (McGraw 
Hill Book Company), 1953. 



12 



TABLES 



Table 1. Reduced matrix elements of quadrupole operators in r— and p— spaces. 
Note that (M \\q p \\ Nl) = (N£' \\q r \\ N£) , ((N ± 2) £' \\q p \\ Nl) — — ((N ± 2) £' \\q r \\ N£). 

(Nl' \\q p \\ Nl) ((N - 2) £' \\q p \\ N£) ((N + 2) £' \\q p \\ Nl) 



pi - p _( N >A Ml+l)(2l+l) / n(n+l+l/2)l(l+l)(2l+l) / 

t ~ 1 ' V 2 / V (2^-l)(2^+3) Y (2^-1) (2^+3) V 



(n+l)(n+l+3/2)l(l+l)(2l+l) 
(2£-l)(2^+3) 



' 6 (n+ 


•l)(n+M-l/ 


'2)£(£- 


-1) 




2£-l 






' 6 n(n 


+£+3/2)(£- 


H)(£- 


H2) 



nl_n_cy / fi (n+l)(n+l+l/2)l(l-l) / 3 (n+l+l/2)(n+l-l/2)l(l-l) / 3 (n+l)(n+2)£(£-l) 

w u ■" ; V 2 2t-\ V 2 2^-1 

/ 3 n(n-l)(l+l)(l+2) /3 (n+l+5/2)(n+l+3/2)(l+l)(l+2) 

V 2 2^+3 V 2 2^+3 
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